The normal contact of rough surfaces with fine discretization

Summary The numerical treatment of a 3-d frictionless normal contact of rough bodies made of a linear–elastic, isotropic material is considered. First, the governing equations and their discretized forms are presented. Due to the order of the problem, the computation time rapidly increases with the grid refinement. Therefore, common computer programs have to be limited to 50×50 contact points. Several numerical techniques (multigrid, iterative solver, distributive relaxation) are used. Cell-centered grid transformation is applied and grid operators have been established for this purpose. In order to handle upto 1000 × 1000 cells the program CON3D has been developed.     

The influence of free convection on soil salinization in arid regions

Evaporation of groundwater in a region with a shallow water table and small natural replenishment causes accumulation of salts near the ground surface. Water in the upper soil layer becomes denser than in the depth. This is a potentially unstable situation which may result in convective currents. When free convection takes place, estimates of the salinity profile, salt precipitation rate, etc., obtained within the framework of a 1-D (vertical) model fail.Very simplified model of the process is proposed, in which the unsaturated zone is represented by a horizontal soil layer at a constant water saturation, and temperature changes are neglected. The purpose of the model is to obtain a rough estimate of the role of natural convection in the salinization process.A linear stability analysis of a uniform vertical flow is given, and the stability limit is determined numerically as a function of evaporation rate, salt concentration in groundwater, and porous medium dispersivity. The loss of stability corresponds to quite realistic Rayleigh numbers. The stability limit depends in nonmonotonic way on the evaporation rate.The developed convective regime was simulated numerically for a 2-D vertical domain, using finite volume element discretization and FAS multigrid solver. The dependence of the average salt concentration in the upper layer on the Rayleigh number was obtained.List of Main Symbols horizontal wavenumber - _L , _T dispersivities (longitudinal and transversal) - D ^* diffusion coefficient (in a porous medium) - g acceleration of gravity - H thickness of the vadoze zone - k permeability - p pressure - Pe Péclet number - q mass flux - Ra Rayleigh number Greek _L , _T dimensionless dispersivities - coefficient of concentration expansion - coefficient of viscosity variation - volumetric fraction of the liquid phase - viscosity - density - stream function - mass fraction of salt in water Vectors and tensors D dispersion coefficient - e unit vector - I unit tensor - J nonadvective salt flux - V liquid phase velocity - x radius-vector     

基于SIMPLER算法的多重网格方法研究

以二维方腔顶盖驱动流为模型,将多重网格方法和SIMPLER算法进行耦合,对不同雷诺数下多重网格加速SIMPLER算法和SIMPLER算法的计算效率进行了对比,数值计算表明:多重网格加速SIMPLER算法不仅能够解决SIMPLER算法不能准确模拟较高雷诺数流场的问题,而且其计算效率远远高于SIMPLER算法.本文也对松弛因子的选取、多重网格实现形式以及网格层数对多重网格加速SIMPLER算法的影响进行了研究,从而为多重网格加速SIMPLER算法的实施提供了计算技术.     

多体系统动力学设计灵敏度分析直接微分法

针对受完整约束的多体系统动力学微分／代数方程数学模型动态最优化设计问题，建立了通用的目标函数和约束方程，并以此为基础，用直接微分方法系统地推导出了计算设计灵敏度的通用公式，最后通过平面机械臂模型对理论结果和相应算法进行了验证．     

用代数应力模型模化型腔充填过程中的紊流流动

针对型腔充填过程中的紊流流动，用代数应力模型研究紊流流动现象．成功地将代数应力模型引入ＰＨＯＥＮＩＣＳ软件中，完成了恒温流体充填过程三维时均速度和自由表面的计算模拟．与实验结果相比说明，本文基于ＰＨＯＥＮＩＣＳ软件开发的充填模拟具备了对复杂型腔充填过程的数值计算．     

An explicit algebraic Reynolds stress and heat flux model for incompressible turbulence: Part I Non-isothermal flow

Tensor representation theory is used to derive an explicit algebraic model that consists of an explicit algebraic stress model (EASM) and an explicit algebraic heat flux model (EAHFM) for two-dimensional (2-D) incompressible non-isothermal turbulent flows. The representation methodology used for the heat flux vector is adapted from that used for the polynomial representation of the Reynolds stress anisotropy tensor. Since the methodology is based on the formation of invariants from either vector or tensor basis sets, it is possible to derive explicit polynomial vector expansions for the heat flux vector. The resulting EAHFM is necessarily coupled with the turbulent velocity field through an EASM for the Reynolds stress anisotropy. An EASM has previously been derived by Jongen and Gatski [10]. Therefore, it is used in conjunction with the derived EAHFM to form the explicit algebraic model for incompressible 2-D flows. This explicit algebraic model is analyzed and compared with previous formulations including its ability to approximate the commonly accepted value for the turbulent Prandtl number. The effect of pressure-scrambling vector model calibration on predictive performance is also assessed. Finally, the explicit algebraic model is validated against a 2-D homogeneous shear flow with a variety of thermal gradients. Dedicated to the memory of the late Professor Charles G. Speziale of Boston University     

求解复杂湍流的非线性涡黏性系数模型和代数应力模型

非线性涡黏性系数模型和代数应力模型联系了线性涡黏性系数湍流模型和完整的微分 雷诺应力模型.随着它们受到日益关注,其形式也越来越多样化.本篇综述的目的是对这些模 型加以总结并比较它们之间的共同点及不同之处,指出它们与完整微分雷诺应力模型之间的 关系,以及相对于线性涡黏性系数模型而言它们在预报流场上所具有的优势.     

The Opposing Effect of Viscous Dissipation Allows for a Parallel Free Convection Boundary-Layer Flow Along a Cold Vertical Flat Plate

External free convection boundary-layer flows are usually treated by neglecting the effect of viscous dissipation. This assumption always results in a non-parallel flow, besides a strong parallel component also a weak transversal component of the (steady) velocity field occurs. The present paper shows, however, that the weak opposing effect of the buoyancy forces due to heat release by viscous dissipation, can give rise along a cold vertical plate adjacent to a fluid saturated porous medium to a strictly parallel steady free convection flow. This boundary-layer flow is described by an algebraically decaying exact analytical solution of the basic balance equations.     

本刊英文版Vol.30(2014),No.9论文摘要(英文)

<正>Least Number of Periodic Points of Self-maps of Lie Groups Jerzy JEZIERSKI Abstract There are two algebraic lower bounds of the number of n-periodic points of a self-map f:M→M of a compact smooth manifold of dimension at least 3:NF_n(f)=min{#Fix(g~n);g~f;g is continuous}and NJD_n(f)=min{#Fix(g~n);g~f;g is smooth}.In general,NJD_n(f)may be much greater than NF_n(f).If M is a torus,then the invariants are